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ABSTRACT 


Our  title  describes  a  phenomenon  best  illustrated  by  our  Theorem  1. 

1.  Let  h(x)  be  an  entire  function  of  exponential  type  A  <  2*.  We 
show  that  there  is  a  unique  entire  function  f(x)  of  exponential  type  A 
satisfying  the  functional  equation 

(1)  f(x+1)  -  f(x)  ■  h(x),  with  f (0)  «  0  . 

2.  We  define  a  function  Sn(x)  which  reduces  to  a  polynomial  in  [0,1) 

(2)  S  (x)  ■  a,x  +  «*_x  +. ..+  a  x  if  0  <  x  <  1  , 

n  1  2  n  ® 

and  we  extend  the  definition  of  Sn(x)  to  all  real  x  by  ashing  that  it 
satisfy 


(3)  Sn(x+1)  -  Sn(x)  *  h(x)  for  all  real  x. 
If  we  require  also  that 

(4)  S(V)(+0)  =  S(V)(-0)  for  v  *  0,1,...,n-1 

n  n 

then  we  have 

(5)  S  (x)  e  Cn_1(R) 

n 

and  also 


(6)  lim  S  (x)  *  f(x)  uniformly  for  all  real  x  . 

n~  n 
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our  tin.  stresses  the  phenceuon  thet  the  c"-'-«,„ti„„ity  rec,„ire„e„t 
(5)  implies  the  approximation  property  (6). 

As  a  second  example  we  deal  with  the  functional  epuation 

(?>  f *  xf(x+1),  with  f(i)  »  i 

satisfied  by  f(x)  *  i/r(x). 
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SIGNIFICANCE  and  explanation 


Examples  are  qiven  of  functional  equations,  like  f(x+1)  -  f(x)  =  h(x), 
where  a  solution  f(x)  is  everywhere  defined,  provided  that  f(x)  is 
prescribed  in  the  interval  0  <  x  <  1.  We  define  a  function  Sn(x)  by 
settinq  Sn(x)  =  Pn(x)  in  0  <  x  <  1,  where  Pn(x)  is  an  as  yet  unknown 
polynomial  of  deqree  n,  and  requirinq  Sn(x)  to  satisfy  our  functional 
equation  for  all  real  x.  We  now  impose  the 

Conti nuity^Requirement.  We  require  the  composite  function  Sn(x)  to 
have  n  -  1  continuous  derivatives  at  the  point  x  *  0. 

It  is  shown  in  our  examples  that  the  Continuity  Requirement  has  the 
followinq  consequences* 

1°,  Sn(x)  has  at  least  n  -  2  continuous  derivatives  for  all  real  x. 

2°.  As  n  ♦  •  Sn(x)  converqes  to  a  solution  f(x)  of  our  functional 
equation. 

The  fact  that  the  Continuity  Requirement  implies  the  approximation 
property  2°  explains  the  meaninq  of  our  title. 


The  responsibility  for  the  wordinq  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


HIGH  ORDER  CONTINUITY  IMPLIES  GOOD  APPROXIMATIONS 
TO  SOLUTIONS  OF  CERTAIN  FUNCTIONAL  EQUATIONS 


T.  N.  T.  Goodman,  I.  J.  Schoenberg  and  A.  Sharma 
1.  Introduction  and  main  results*  We  describe  here  an  apparently  new 
principle  of  generating  polynomial  approximations  in  0  <  x  <  1  to  the 


solutions  of  certain  functional  equations* 

Our  mysterious  title  will  become  less  so  if  we  mention  first  the  example 
of  the  exponential  Euler  spline  Sa(x;t)  satisfying  the  functional  equation 

(1.1)  Sn(x+1;t)  =  tSn(x;t),  with  Sn<0;t)  «  1  . 

Here  t  is  a  constant  such  that 

(1.2)  t  =  1 1 |e10t,  -w  <  a  <  w,  t  +  0,  t  j*  1  . 

Let 

(1.3)  pn(x)  »  1  +  CjX  +  CgX2  +...+  cnxn 

be  an  as  yet  unspecified  polynomial.  We  define  the  function  SR(x;t)  by 
requiring  firstly,  that 


(1.4)  S  (x;t)  =  P  (x)  if  0  <  x  <  1  , 

n  n  * 

secondly,  that  Sn(x;t)  should  satisfy  the  functional  equation  (1.1)  for  all 
real  x.  The  result  so  far  is  that  Sn(x;t)  is  a  piecewise  polynomial 
function  depending  on  the  n  parameters  c^,...,c  . 

Now  comes  the  essential  requirement:  We  ask  that 

(1.5)  Sn(x;t)  e  Cn-1(R)  . 

This  is  what  turns  Sn(x;t)  into  the  exponential  Euler  spline.  In  fact  it 
suffices  to  require  Cn_1 -continuity  only  near  the  point  x  =«  0,  because 
(1.1)  will  propagate  this  order  of  continuity  to  all  other  integer  points. 
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Concerning  the  continuity  at  x  =  0:  By  (1.4)  and  (1.1)  we  see  that  in 


0  <  x  <  1  we  have  Sn(x;t)  *  Sn(x)  *  pn(x),  while  in  -1  <  x  <  0  we  have 
Sn(x)  =  Sn(x+1)/t  =  pn<x+1 )/t.  These  two  functions  Pn(x)  and  pn(x+1)/t 
will  join  at  x  *  0  with  n  -  1  continuous  derivatives  iff 


(v)  (v) 

Pl  '(0)  =  P*  (1)/t  or 
n  n 


(1.6) 


P<V)(1)  *  tP<V>(0),  (v  =  0, . . . ,n-1 )  . 

n  n 


These  are  precisely  the  equations  that  define  the  polynomial 


(1.7) 


pn(x)  =  An(x;t)/An(0;t)  , 


where  An(x;t)  is  defined  by  Euler's  generating  function 


(1.8) 


„  „  “A  ( X ; t ) 

t-1  xz  r  n  n 

- e  -  l  - - -  z 

.  z  «  nl 

t-e  0 


For  details  see  [9],  where  the  exponential  Euler  splines  were  first 


introduced.  It  was  also  shown  there  that 


(1.9) 


lim  S  (x;t)  =  t*  *  It^e*0*  for  all  real  x  . 
n-*»  " 


To  summarize:  We  assume  (1.1)  and  (1.4);  now  the  continuity  condition 
(1,5)  implies  the  approximation  property  expressed  by  the  limit  relation 


(1.9). 


That  continuity  under  similar  circumstances  implies  approximation  is  not 
an  isolated  phenomenon.  Our  second  example  is  the  functional  equation 


(1.10) 


f(x+1)  -  f(x)  *  h(x),  with  f(0)  *  0  , 


satisfied  by  the  "sum"  f(x)  of  a  prescribed  function  h(x).  We  assume 
h(x)  to  be  entire  of  exponential  type. 

Our  main  result  concerning  (1.10)  is 
Theorem  1.  1.  We  assume  that 

nasasasa  —  ■  ■"  —  ■—  ■■  1 

.  (v)  i  v 

(1.11)|h  ( 0 ) |  <  C» A  ,  (v  -  0,1,...  ),  where  C  and  A  are  constants 

such  that 


A  <  2»  . 


(1.12) 


The  equation  (1.10)  has  a  unique  entire  solution  f(x)  such  tha t 
(1.13)  |f(V>(0)|  <  K*AV,  (v  »  0,1,...  ),  with  K  constant  . 

2.  We  define  Sn(x)  by  setting 

d.14)  v*>  mTTai*  ♦srv2  ♦•••■ *sr v"  0 i* <  •  • 

and  extend  its  definition  to  all  real  x  so  as  to  satisfy 

(1.15)  Sn(x+1)  -  Sn(x)  -  h (x )  . 

If  we  require  that 

(1.16)  Sn(x)  is  of  class  C*1”1  in  a  neighborhood  of  x  =  0  , 

then 

(1.17)  Sn(x)  €  c"-1  (R) 

Finally 

(1.18)  lirn  S  (x)  =  f(x)  uniformly  on  R  . 

n~  n  - 

For  a  proof  see  42.  Theorem  1  is  very  close  to  a  result  of  J.  M. 
Whittaker  [12,  Theorem  3  on  paqe  22],  Our  proof  hrinqs  out  more  clearly  the 
main  idea  of  this  paper:  That  (1.14),  (1.15),  and  (1.17)  imply  (1.18). 

If  e.q.  h(x)  =  231  (satisfyinq  (1.11)  with  A  *  loq  2  <  2*),  then  the 
unique  solution  of  (1.10)  satisfyinq  (1.13)  is,  of  course,  f(x)  *  2X.  The 
approximation  Sn(x)  to  2*,  furnished  by  Theorem  1,  is  different  from  the 
exponential  Euler  spline  Sn(x;t),  for  t  =  2. 

Our  third  and  last  example  is  perhaps  the  most  interesting  one:  We 
consider  the  functional  equation 

(1.19)  f(x)  *  xf (x+1 ) ,  with  f(1)  *  1  , 

satisfied  by  the  reciprocal  T-function 


Notice  that  (1.19)  implies  f'(x)  *  f(x+1)  +  xf'(x+1),  and  settinq  x  =  0  we 
obtain  that 


(1.21) 


f  *  (0)  =  f  ( 1 )  *  1 


Also  (1.19)  implies  for  x  *  0  that  f(0)  ■  0.  Accordingly,  we  consider  the 


polynomial 


(1.22) 


P_(x)  «  x  +  a,x2  +...+  a_xn  , 


and  following  our  general  approach  we  define  the  function  sn(x)  by  setting 


(1.23) 


S  (x)  *  P  (x)  if  0  <  x  <  1  , 

n  n  * 


and  extend  its  definition  so  as  to  satisfy 


(1.24) 


Sn(x)  *  xSn(x+1),  sn(1)  *  1,  for  all  real  x 


We  pass  now  to  the  critical  continuity  reguirements  on  Sn(x).  Observe 
that  by  (1.23)  and  (1.24)  we  obtain  the  following: 


(1.25) 


In  0  <  x  <  1  we  have  S  (x)  ■  P  (x)  , 

*  n  n 


(1.26)  in  -1  <  x  <  0  we  have  S  (x)  *  xS  (x+1 )  «  xP  (x+1)  . 

»  n  n  n 

If  we  want  these  two  functions  to  join  at  x  *  0  with  n  -  1  continuous 
derivatives,  we  must  ask  that 


(1.27) 


P  (x)  -  xP  (x+1)  *  0(x  ),  as  x  +  0  . 

n  n 


Our  main  result  is 


Theorem  2.  For  n  =  1,2,...  there  is  a  unique  polynomial 


(1.28) 


_  ,  v  (n)  2  (n)  n 

P  (x)  *  x  +  a „  x  +...+  a  x 
n  2  n 


satisfying  (1.27)  with  Pn (1 )  *  1 .  Its  coefficients  are  rational  numbers. 

A  proof  is  given  in  (3.  The  first  six  polynomials  of  this  remarkable 

sequence  *re 

Pi (x)  *  Pj(x)  *  x 

_  ,  .  A  1  2  13 

P3(x)  *  X  +  —  x  -  —  X 

(1.29)  *  x  +  —  x  ~  f  *  +  J  * 

_  ,  .  42  23  31  425 

P5(X)  -X  *jx  -^x  ♦  35  *  +  35* 


152  173  1  4  4  5  1  6 

Pfi(x)  *  x  +  26  X  26  X  26  X  +  26  X  26  X 
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The  question  as  to  how  the  Cn“^ -continuity  of  Sn(x)  at  x  =  0  is 
transmitted  by  (1.24)  to  the  other  integers  is  easily  answered  as  follows. 

(1.30)  In  f-2,-1)  Sn(x)  =xSn(x+1)  =  x(x+1 )Sn(x+2)  =  x(x+1 )Pn(x+2)  . 

From  (1.27)  we  find  that  Pn(x+1)  -  (x+1)Pn(x+2)  =  0(x+1 )n  as  x  ♦  -1, 
whence 

xP  (x+1 )  -  x(x+1)P  (x+2)  =  0x(x+1)n  =  0(x+1)n  as  x  ♦  -1 
n  n 

But  then  (1*26)  and  (1*30)  show  that 

S(V)(-1+0)  *  s(V)(-1-0),  (v  =  0, . . . ,n-1 )  , 

m  n 

and  we  readily  find  similarly  that 

( 1 . 31 ) S^ (k+0)  =  S^V^(k-0)f  (v  s  0, ...,n-1)  for  all  integers  k  <  0  . 

n  n  ■■  1  ■  — — 

Matters  are  slightly  different  if  we  pass  to  positive  k. 

S  (x-1 )  P  (x-1 ) 

(1.32)  In  M,2)  :  S  (x)  =  — - ; -  =  — - ; -  . 

n  x-1  x-1 

Now  (1.27)  shows  that  Pn(x-1)  “  (x-1)Pn(x)  =  0(x-1  )n  and  so 

P  (x-1) 

(1.33)  — —  - Pn^x^  *  0(x-1)n  as  x  1  . 

But  then  (1.23)  and  (1.32)  show  that 

(v)  (v) 

S'  '(1+0)  =  S'  '(1-0)  only  for  v  =  0,...,n-2  , 

n  n  - 

so  that  we  have  lost  at  x  =  1  one  order  of  continuity. 


From  here  on  there  is  no  further  loss  of  continuity: 


(1.34) 


In  [2,3)  :  S  (x) 
n 


1  S(x-2)  =  _1_  P(x-2) 
x-1  x-2  “  x-1  x-2 


Pn(x_2)  -1 

However,  (1.33)  gives  -  Pn(x-1)  =  0(x-2)n  as  x  *  2,  whence 


This  continues  for  all  positive  k  and  we  obtain 


Corollary  1.  The  approximation  S_(x)  defined  by  (1.23),  (1.24)  and 

SSS333S5S  1  ■  ■  '  '  1  '■  ■  '  1  *1  "  ' 

(1.27)  has  the  following  continuity  properties: 

(1.35)  sn(x )  1  Cn_1  at  the  points  x  *  0,-1, -2,...  , 

(1.36)  Sn(x)  e  c””2  at  the  points  x  =*  1,2,3,... 

The  only  exception  is  S^x)  which  is  continuous  on  all  of  R. 

The  exceptional  case  of  S^x)  is  due  to  P.,(x)  =  x  and  the  identity 
x  -  x(x+1)  =  x^  =  Ox ^  as  x  +  0. 

A  further  immediate  consequence  of  (1.24)  is 

Corollary  2.  The  approximation  Sn(x)  interpolates  1/T(x)  at  all 
integers,  hence 

(1.37)  Sn(V)  *  r|v)  f°r  V€  Z  • 

We  have  referred  above  to  Sn(x)  as  an  approximation  of  1/r(x).  Is  it? 

This  essential  point  we  can  not  settle  and  wish  to  state  it  as 

Coniecture  1.  The  sequence  of  polynomials  P„(x)  of  Theorem  2  satisfy 

3=3333=333  1  ■  — —  '  1  ■l,,‘  1»  "  ' 

(1.38)  lim  P  (x)  =  n uniformly  in  every  circle  of  <C  . 

n-n»  n  r(x) - - 

As  an  extremely  weak  substitute  we  establish 
Theorem  3.  For  n  »  1,2,...  we  have 

(1.39)  P  (x)  >  0  if  0  <  x  <  1  , 

n  —  = 

(1.40)  P'(x)  >0  if  0  <  x  <  1  . 

n  «  = 

Proof  of  (1.39);  We  apply  the  Budan-Fourier  theorem  in  its  classical 
formulation;  see  e.g.  (8)  where  it  is  derived  from  the  Descartes- Jacobi 
theorem  by  means  of  a  special  totally  positive  matrix;  see  also  [5,  316-318]. 
Writing  Pn(x)  ”  F(x)  we  use  the  notation 

(1.41)  V(x)  -  v(P(x),P’(x),...,P(n)(x)) 

for  the  nu  ber  of  inges  of  sign  in  the  sequence  as  indicated.  If  Z(0, 1)  * 

number  of  zetos  of  P(x)  in  (0,1),  then  the  theorem  states  that 

(1.42)  Z(0,1)  <  V(1)  -  V(0)  . 


-6 


I 


However,  the  relation  (1.27)  will  provide  much  information  concerning  V(0) 

and  V(1):  Writinq  F(x)  =  P(x)  -  xP(x+1)  we  find  by  Leibnitz's  formula  that 

(v)  (v)  (v)  (v-1 ) 

F  (x)  =  p'  '(x)  -  xP'  ( x+1 )  -  VPV  (x+1).  Settinq  x  =  0,  (1.27)  shows 


that 

(1.43) 


(v)  (v-i) 

P  '(0)  =  vpv  '(D,  (v  =  1,2,...,  n-1 )  , 


and. the  riqht  side  of  (1.42)  promises  to  vanish.  In  fact  from  P(x)  = 

x  +  a^x^  +. ..+  anxn  and  (1.43)  we  find  that 

0.44)  V(0)  =  v(  0, 1  ,a2, . . .  ,an_2,  an_i»  aR) 

(1.45)  V(1)  =  v(1,a2,a3,...,an_1,  P(n_1)(1),  an)  . 

We  distinquish  two  cases. 

1.  an_ian  <  0.  In  this  case  we  see  from  (1.45)  that  V(1)  =  V(0), 

whatever  the  siqn  of  P^n-1^(1)  miqht  be. 

2.  a„  >  0.  In  this  case 

n—  I  n 

p(n-1)(x)  _  an_1(n-1)l  +  ann!x 

shows  that  P^n-1^(1)  has  the  common  siqn  of  and  an,  so  that  aqain 

V( 1 )  =  V( 0)  . 

Above  we  have  assumed  that  an_ian  7*  0,  but  the  last  result  is  seen  to 
hold  in  any  case,  whether  one  or  both  an_i  and  an  should  vanish. 

We  postpone  a  proof  of  (1.40)  to  (4  because  it  used  results  from  $  3. 
Usinq  (1.24)  as  in  our  proof  of  Corollary  1  we  obtain 
Corollary  3.  We  have 

(1.46) 
and 

(1.47)  sn(x)  <  0  lH  (-1.0),  Sn(x)  >  0  in 

(-2,-1),  Sn(x)  <  0  in  (-3,-2)  a.s.f. 

with  alternatinq  siqns. 


sn(x)  >  0  II  x  >  0  , 


A  curious  neqative  consequence  is 


-7- 


Corollary  4.  The  function  G(x)  ■  loq  Sn(x)  is  not  a  concave  function 
for  x  >  0. 

This  follows  from  a  theorem  of  Bohr,  Mollerup  and  Artin  f2,  Theorem  2.1 
on  paqe  141  that  F(x)  =  loq  T(x)  is  the  only  convex  solution  of  the 
functional  equation 

F(x+1)  -  F(x)  =  loq  x  for  x  >  0,  with  K ( 1 )  =  0  . 

It  suffices  to  observe  that,  by  (1.24)  and  (1.46),  the  function  G(x)  = 
loq  Sn(x)  satisfies 

G(x+1 )  -  G(x)  =  -loq  x  for  x  >  0,  with  G(1)  »  0  , 

while  the  identity  G(x)  =  -F(x)  is  evidently  impossible. 

Everybody  believes  the  Euler  constant 

n 

Y  =  lim(I  —  -  loq  n)  =  .57721  56649  ... 

1  v 

to  be  irrational,  but  as  far  as  we  know  nobody  has  proved  this.  We  also 
believe  it,  and  will  here  assume  that 

(1.48)  Y  is  an  Irrational  number  . 

By  (1.37)  we  know  that  SR(x)  interpolates  T  \x)  at  all  integers. 
However,  much  more  is  true  as  stated  in 

Corollary  5.  We  assume  (1.4R)  to  hold,  and  let  n  >  3.  For  each  k  e  Z 

3Ba  =  3*2I33SSi3:  - -  1  "  1  ■  ■  *  1,1 

the  two  open  arcs 

(1.49)  V  =  S  (x)  (k  <  x  <  k+1 )  and  y  *  T  1 (x)  (k  <  x  <  k+1 ) 

n  ■ 

cross  each  other. 

Proof.  Let  k  *  0.  Assuminq  that  the  arcs  (1.49)  do  not  cross  and  let 
us  qet  a  contradiction.  More  precisely  let  us  assume  that 

(1.50)  S  (x)  -  P  (x)  <  r_1(x)  if  0  <  x  <  1  . 

n  n  « 

From  (1.32)  we  have 


S  (x) 


(x— 1 ) 

■j - #  r~  (x) 


r-1(x-i ) 

X— 1 


if 


1  <  x  <  2  , 


-8- 


I 


it 

and  therefore 

S  ( x '  S  (x-1) 

- -  =  -  <1  by  (1.50)  . 

r_1(x)  r-1(x-1)  ' 

Hence 

(1.51)  s  (x)  <  r‘\x)  if  1  <  x  <  2  . 

n  = 

However,  (1.50)  and  (1.51)  show  that  S  (x)  <  T  ' (x )  throughout  the  interval 

n  — 

0  <  x  <  2*  Since  S  (1)  =  T  ^  ( 1  >  #  while  S  c  c\r),  we  conclude  that  their 

n  n 

slopes  at  x  =  1  must  be  equal*  This  is  impossible  by  (1.48)  because 

(d/dx)r_1(i)  =  -r*(1)/r(1)2  =  -r ' ( I )  =  Y,  while  S'(1)  is  evidently  a 

n 

rational  number. 

We  illustrate  numerically  our  conjecture  (1.38)  as  well  as  Corollary  5 

for  the  case  when  n  =  6.  With  Pg(x)  as  given  in  (1.29)  and  the  table  of 

T(x)  f1.  Table  6.1  on  page  267)  we  find  that  the  functions  Pg(x)  and 

T  \x)  cross  at  a  point  €  ■  .1182...,  and  that  we  have 

0  <  r_1(x)  -  P,(x)  <  .000  001  if  0  <  x  <  C  , 
o 

0  <  P  (x)  -  r_,(x)  <  .000  15  if  C  <  x  <  1  . 

o 

From  the  evident  proportion 

S  (x)  S  (xfk) 

=  «  •  (x  ^  2,  K  €  X)  , 

r  (x)  r~  (x+r) 

it  is  clear,  that  for  n  =  6,  the  two  arcs  (1.49)  cross  at  the  point 
x  =  . 

Concluding  we  wish  to  point  out  that  the  two  papers  (4)  and  [11 J  deal 
with  subjects  related  to  the  topic  of  the  present  note. 
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2.  The  equation  f (x+1 )  -  f(x)  =  h(x):  ^_proof _of _Theorem  1.  We  apply 
our  previous  method:  let  us  construct  a  compos i tie  function  Sn(x)  having 
the  two  properties: 

(2.1)  S(x)=P(x)=-7-a*x+4ra~x2+.  ,.+  ^— axn  if  0<x<1  , 

n  n  1!  1  2!  2  n!  n  = 

(2.2)  S  (x)  satisfies  S  (x+1)  -  S  (x)  =  h(x)  for  all  real  x  . 

n  n  n 

On  this  function  we  now  impose  the  continuity  condition: 

(2.3)  S  (x)  e  Cn_1(R)  . 

n 

We  claim:  _If_  Sn(x)  is  of  class  C0-1  in  a  neighborhood  of  x  =  0 
then  (2.3)  holds. 

Proof.  We  write  S  (x)  =  S(x).  Differentiation  of  (2.2)  shows  that 
- — —  n 

(V)  (V)  (V) 

Sl  '  (k+1+0)  -  s'  '  (k+0)  =  h'  '(k+0), 

S{V)(k+1-0)  -  S(V)(k-0)  =  h(V)(k-0) 

and  by  subtraction 

(2.4)  S(V,(k+1+0)  -  S(V)(k+1-0)  =  S(V,(k+0)  -  S(V)(k-0)  . 

This  implies  that  if  the  riqht  side  of  (2.4)  vanishes  for  k  «  0,  it  must 
vanish  for  all  k  and  (2.3)  is  established. 

To  inforce  Cn_1 -continuity  of  Sn  at  x  *  0  we  observe  the  followinq: 

By  (2.1)  and  (2.2)  we  have: 

In  (0,1)  :  S(x)  =  Pn(x),  while  in  (-1,0)  :  S(x)  =  S(x+1)  -  h(x)  =  Pn(x+1) 

h(x).  Therefore  the  equations  S^V\+0)  *  S^V\-0),  (v  =  0,...,n-1) 

(v)  (v)  (v) 

are  equivalent  to  P  (0)  »  P  (1)  -  h  (0),  and  hence  to 

n  n 

(v)  (v)  (v) 

(2.5)  Pl  '(1)  -  P'  '(0)  +  h'  ;(0),  (v  .  0,...,n-1)  . 

n  n 

The  remainder  of  the  proof  is  divided  into  several  parts, 
a)  Construction  of  the  Pn(x).  In  terms  of  our  explicit  expresssion 
(2.1)  of  Pn(x),  the  system  (2.5)  is  seen  to  be 
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(2.6) 


For  the  matrix  of  (2.6)  we  find  that 


and  therefore 


1  1  1  n-1  e  -1  _ ,  n . 

t  -  —  +  0<t) 


n-1  B 


.  ..1  1  1  n-1 ,  t  n.  r  v  n. 

i/(tt  +  ir1  +,"t7n-t  >  --rr  +  0<t  >  -  l  vr  *°(t  > 


where  By  are  the  Bernoulli  numbers.  It  follows  that  the  solution  of  (2.6) 
is  qiven  by  the  matrix  equation 


1 

B1 

h(  0) 

a1 

1! 

(n-1)! 

0 

1 

...  n-2 

• 

h*(0) 

• 

a2 

s 

(n-2) ! 

e 

• 

• 

e 

• 

m 

* 

an 

0 

• 

1 

o 

7 

c 

7c 

or 


and  so 


,  nyk  Bv  (k+v-n 


(0) 


n  k  n  k  n-k  B%1  „  ... 

b/n  v  x  r  x  r  v  (k-1+v) 

n  k-1  k!  k  k-1  kl  v-o  V! 


(0) 


Replacinq  k  by  k  +  1  we  find  that 

(2.9) 


n-1  k+1 

P (x)  -  l  X 


n 


k=0 


where 

(2.9*) 


n,k 


"'y'1  Bv  (v-rit) 

v*o 


(k+1)!  n,k 


n-1  B 


(0) 


T  v~*  h(v)fm 
vf^  (v-k)!  h  (0) 


6)  Construction  of  the  solution  f(x).  Let  m  show  that  the 
polynomials  (2.9)  converqe  if  0  <  x  <  1  and  determine  their  limit.  From 
lim  sup  |bv/vi|1/V  *  1/2 x  we  conclude  that  for  a  small  6  >  0  we  have 
|bvi/vj|1/V  <  l/(2»-$)  for  sufficiently  larqe  v,  and  so 

Ib. 


V 

V! 


<  C, ( for  *11  v  • 


Rut  then  our  assumption  (1.11)  shows  that 

|B 


(2.10) 


V  fV+le)  I  k  A  \J 

—  h'  ^  ( 0)  |  <  C0A  (  )  for  all  v,  with  constant  . 


| vi . |  '  2  2w-o  —  r  2 

In  view  of  (1.12)  we  may  assume  6  to  be  so  chosen  that  A  <  2*  -  5.  Then 


(2.10)  shows  that  is  well  defined  by 

*  B 


(2.11) 


I  wh 

K  V-0  1 


V  .  (v+k) 


(0) 
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and  that 


(2.12) 
But  then 
(2.13) 


| S.  |  <  for  all  k,  constant. 


f(x)  =  l 


k+1 


k=0 


(k+1)  l  k 


is  an  entire  function  satisfyinq  ( 1 • 1 3 ) • 

Y)  Let  us  show  that 

(2.14)  lim  P  (x)  =  f(x),  uniformly  in  0  <  x  <  1  • 

n  1  ■  -  s  ss 

n+* 

By  (2.9)  and  (2.13) 

n-1  k+1  08  k+1 

f  (x)  -  Pn(x)  =  I  (k+1 )!  (Sk  "  An,k)  +  (k+1  )!  Sk  ' 

k=0  k«n 

and  that  it  suffices  to  assume  0  <  x  <  1,  and  to  show  that 


n-1  k+1 

R  (x)  »  Z  X 


n'  (k+1 )  1  k  n,k 


uniformly  in  f0,1]. 

In  [0,1]  we  have  from  (2.10)  that 


(S,_  -  A  .  >  -*•  0  as  n  ♦  -  , 


,pn(x)|  -  ^  TkTTTT  5  lyf  h('Mk)(0,l 
n  =  k=0  {k+1)!  v=n-k  V! 


n-1 


is  i 


1 


I 


v+k 


k»0  (k+1)!  v*n-k  (2ir-6) 


and  replacinq  the  index  v  by  v+n-k 

00 

I 


n-1  •  v+n 

,Rn(x)l  *  C2  E  nrrnr  2  - 

2  k=*0  (K+1)  v*0  (2*-$)  K+ 


.  n-1  *  ...  r.k  _ 

„  „  A  xn  v  V  (2W-0)  ,  A 

J„  i0  <2^' 


A  ,n  v  V  (2n)  .  A  .v  .  _  ,  A  .n 


<  Co(^TT)n  l  l 


(^7)  <  CA 


4'2i^I 


which  establishes  (2.14). 

6)  Proof  that 

(2.15)  f(x+1)  -  f(x)  *  h(x )  . 

Let 

(2.16)  CL(x)  *  Pn(x+1)  -  P„(x)  . 

n  n  n 

(v)  (v)  (v)  (v) 

By  (2.5)  we  have  O  '  (0)  =  P  (l)-P'  7 (0)  =  h  (0),  which  shows  that 

n  n  n 


Q  (x) 

n 


T  ^h<v,(o) 

o 


But  then  (2.16)  becomes 

"v1  xV  (V) 

P  (x+1 )  -  P  (x)  *  I  — t-  h '  '(0) 
n  n  Q  vi 

and  now  (2.14)  implies  (2.15). 

Remarks.  1.  Is  f(x)  the  unique  solution  of  (1.10)  satisfyinq 

(1.13)?  That  it  is  we  see  as  follows:  The  difference  of  two  solutions  of 

(1.10)  is  also  entire  of  exponential  type  <  2*  and  is  also  periodic  of 

period  1.  By  a  qeneral  theorem  (see  [3,  Theorem  6.10.1  on  paqe  109])  this 

difference  must  reduce  to  a  constant. 

2.  Since  S  (x)  and  f(x)  satisfy  (1.15)  and  (1.10)  we  have 
n 

Sn(x+1)  -  Sn(x)  -  (f(x+1)  -  f(x))  =  h(x)  -  h(x)  =  0 

and  so 

S  (x+1)  -  f(x+1)  *  S  (x)  -  f(x)  for  all  x  . 
n  n 

This  shows  that  the  uniform  converqence  in  (2.14)  implies  that 

(2.17)  lim  S  (x)  *  f(x)  uniformly  for  all  real  x  . 

n  — ^ — 
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3*  The  equation  f(x)  =  xf(x+1):  A  proof  of  Theorem  2*  The  matter 

a3S33S333=22  ========= 3333S33S3 

seems  quite  simple  if  we  use  the  proper  tools.  From  (1.27)  for  x  *  0  we  get 

that  P  (0)  =  0.  For  P(x)  =  x  the  left  side  of  (1.27)  becomes  x  -  x(x+1)  = 
n 

-x2,  and  this  shows  that  P-|(x)  =  P2(x)  =  x. 


Assuminq  n  >  2,  (1.27)  qives  P'(x)  -  P  (x+1)  -  xP'(x+1)  =  Onx  ,  and 

n  n  n 

for  x  =  0  that  P'(0)  =  P  (1)  =  1.  For 

n  n 


(3.1 ) 


P  (x)  =  x  +  a,x2  +. ..+  a„xn 
n  2  n 


we  write  (1*27)  explicitly  and  find 

2  3  n 

x+a^x  +  a  -x  +.  •  •  +  a  x 
2  3  n 

-  x(x+1 ) 


-  a^xfx  +  2x+1 ) 

3  2 

-  a^x(x  +  3x  +3x  ♦  1) 


,  n  ,n.  n- 1  ,n.  n— 2  , .  _  n  .  . 

-  a  x(x  +  (,)x  +  (_)x  +...+  1 )  •  Ox  as  x  ♦  0  . 

n  1  2 

Collectinq  terms  and  writinq  that  the  coefficients  of  x,x2, . . . ,xn-1 ,  vanish, 
we  qet  the  equations 

-a2  *  a3  . an-2  -  an-1  *  an  *  0 


a2  -Ja2'3a3 . 'n;2»°„-2  *  <V>V,  -  <>„  *  ° 


_  ,n-2.  ,n-1 .  ,n.  n 

l2  ‘  ^3  - (  2  ,an-2  "  (  2  )an-1  “  (2)an  “  ° 

-  a  . ("“2>a  o  "  1  "  (3)an  “  ° 

3  3  n-2  3  n-i  3  n 


-a  -  (n“]!)a  -  (  ")a  *  0  . 

n-2  n-2  n-1  n-2  n 


This  is  a  system  of  n-1  equations  in  a2'<  Rearranging  these 


equations  we  may  write 
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*  0 


It  is  obtained  from  a  solid  minor 

111  1 
2  3  4  5 

1  3  6  10 

0  1  4  10 

of  the  Pascal  triangle  by  subtracting  1  from  the  elements  of  the  diagonal 
just  below  its  main  diagonal.  A  simple  induction  shows  that  all  these  minors 
are  »  Is 


1  6 


1  1 


2  3 


(3.3) 


Dn  = 


1  3 


1  -  (n:2) 


Indeed  by  successive#  nantractinq  each  column  from  the  next  one  we  find 


0  0 


1  1 


3  4 


=  D.  a.s.f. 
4 


3  6 


We  now  consider  the  determinant  of  the  system  (3.2) 


1  1 


1  3-1 


<";2> 


(3.4) 


C>-  <n-2> 


and  wish  to  prove 


Lemma  1#  We  have 


(3# 5)  A  >  1  for  n  *  1,2? #e#  • 

n  *  — 

Proof.  1°.  We  use  the  fact  first  pointed  out  in  [71  (see  also  r 8 1 )  that 
the  infinite  Pascal  trianqle  is  totally  positive,  i.e.  all  its  minors,  of  all 


orders,  are  >  0. 
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2°*  We  split  the  first  n  -  2  columns  of  (3*4)  into  two  columns,  the 

second  columns  containinq  only  vanishinq  elements,  except  the  sinqle  element  * 

-1.  In  this  way  A  as  a  sum  of  2n~^  determinants 

n 

(3.6)  A  -  D  +  l  A(i  ,i  ....  ,i  )  ,  (1  <  i  <  i  <...<  i  <  n-2)  . 

n  n  r>1  1  2  r  =12  r  = 


where  r>n  =  1,  by  (3.3),  while  (ijf...,i  )  runs  throuqh  all  combination 
of  r  amonq  the  numbers  1,2,...,n-2,  (r>1).  Thus  A(i^,...,i^)  is  the 
determinant  obtained  from  Dn  by  replacinq  its  columns  i.j  *  * « ,ir  by  the 
columns 


(0, 


u.+n 

rO,  -  1  ,  0, 


.0) 


(0, 


(i +1) 

,0,  -ri  ,  0, 


r0) 


We  summarize  its  structure  by 


(3.7) 


Mi 


I'** 


,ir>  = 


(ix)  (i2)  (ir) 


Here  all  elements  are  the  old  elements  of  nn  except  for  those  in  the  r 
columns  ij .ij# • • • »ir« 

3°.  We  apply  Laplace's  rule  (see  (5,  paqe  61)  of  expanding  the 
determinant  (3.7)  by  all  its  minors  of  order  r,  from  its  r  columns 
i^,...,ir,  multiplied  by  their  alqebraic  compliments.  Hi  is  expansion  reduces 


-1  R- 


evidently  to  a  sinqle  term 


-1 


-1 


0 


(3.8) 

A(i j , • • • , i  )  = 

•  i 

• 

• 

• 

c 

X  C  , 

rxr 

where 

r 

r 

(3.9) 

l 

c  =  (-I)n=1 

a  +D  +  l 
i-i 

i  . 

3 

x  P  , 

where 

V  is  a  minor  of  D  -  and 

n 

therefore  V 

>0  by  1°. 

4°.  From  (3*8)  and  (3.9)  we  find  that 


2li  *r 

(3.10)  A(i  f...#i  )  =  (-1)  (-1)  1  •  P  =  P  >  0  . 

Now  (3.3),  (3.6)  and  (3.10)  show  that  (3.5)  is  established. 

Remark.  Elimination  of  the  unknowns  aj, • • • ,an  between  (3.1)  and  the 
system  (3.2)  allows  us  to  write  PR(x)  explicitly  in  terms  of  a  quotient  of 
two  determinants,  e.q. 

2  3  4  5  , 
xxxx  1111 

1111  2-1345 

P_(x)  =  x  +  . 

1  3-1  6  10  *  1  3-1  6  10 

0  1  4-1  10  0  1  4-1  10 

which  is  easily  seen  to  aqree  with  its  expression  as  qiven  in  (1.29).  The 
above  expression  of  Pj(x)  is  readily  continued  for  n  >  5  by  simply  adding 
obvious  last  rows  and  columns  to  both  determinants;  the  only  trouble  is  that 
we  obtain  Pn(x)  =  x  +  a  ratio  of  two  determinants  of  order  n  -  1. 
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4.  A  proof  of  the  second  part  (1.40)  of  Theorem  3.  We  know  that  the 

====33=3=====  333=333=33333333333333332333=3333 

coefficients  aj , ...,an  of  the  polynomial  (3.1)  are  the  solutions  of  the 
system  of  equations  (3.2).  Let  us  use  this  fact  to  show  first  that 


(4.1) 


By  Cramer's  rule 


(4.2) 


*2  -ip;<o>  > 0  • 


S  =  VAn  ' 

2  n  n 


where  A  is  qiven  by  (3.4),  and  A  is  obtained  from  A  by  omittinq  its 
n  n  n 

first  column  and  second  row.  We  know  by  (3.5)  that  A  >1.  Let  us  show  that 

n  —  ■  ■ 

also 


(4.3) 


A  >  1  . 

n  * 


This  is  done  by  applyinq  to  the  determinant  A  precisely  the  procedure 

n 

previously  applied  to  A  to  prove  (3.5).  Here  we  need  the  followinq 

n 

stronqer  form  of  the  total  positivity  of  the  matrix  of  the  determinant  Dn  of 

(3.3): 


( 4. 4 )  A  minor  of  Dn  which  doe  not  vanish  formally,  i.e. 
because  it  has  too  many  zero  elements,  is  positive. 

Denote  by  the  minor  of  Dn  obtained  by  omittinq  its  first  column 

and  second  row.  This  minor  havinq  an  inteqer  value,  the  property  (4.4) 
implies  that 

(4.5)  D  >  1  . 

n  * 

Prom  this  point  on  the  proof  of  (3.5)  as  qiven  in  43  applied  without  any 

essential  chanqe  to  establish  (4.3).  Now  (3.5)  and  (4.2)  imply  (4.1). 

To  establish  (1.40)  we  aqain  apply  the  Budan-Fourier  theorem,  this  time 

to  P'(x).  Let  2'  denote  the  number  of  zeros  of  P'(x)  in  0  <  x  <  1. 
n  n 

Writinq  V'(x)  »  v(P'(x),  P"(x),... ,P^n^(x) )  we  obtain  by  that  theorem  that 

n  n  n 

(4.6)  2*  <  V'(1)  -  V'(0)  , 


while  the  equations  (1.43)  show  that 


V 1 ( 0 )  =  v(1,  a  ,  a  ,  a  ) 

2  n— 2  n-1  n 

v'(1)  =  v(a_,a  ,...,a  P(n_1)(1),  a  )  . 

2  3  n— I  n  n 

Since  a2  >  0  by  (4.1),  it  follows  as  in  our  proof  of  (1.39),  that  the  riqht 

side  of  (4.6)  vanishes,  hence  7,'  =  0.  Since  P'(0)  =  1,  we  have  established 

n 

(1.40). 
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5.  On_tHe  _na  ture  _of  _Con  _of  _(1»382»  The  coefficients 

of  the  expansion 

ao 

(5.1)  P(X)  =  ^  ckx  '  (c1  =  D 

are  qiven  in  M,  paqe  2561  to  15  decimal  places,  as  oriqinally  computed  by 
H.  T.  Davis  in  1933,  with  corrections  due  to  H.  E.  Salzer.  If  we  substitute 
the  expansion  (5.1)  of  f(x)  =  1/T(x)  into  the  functional  equation  f(x)  = 
xf(x-tl)  we  find  that  the  satisfy  an  infinite  system  of  linear  equations 

obtained  from  our  system  (3.2)  by  substitutinq  Cy  for  and  then  lettinq 

n  ♦  °°. 

The  problem  is  to  show  that  the  solutions  (k  -  2,...,n)  of  the 

partial  system  (3.2)  converqe  element-wise  to  the  solutions  of  the 

infinite  system.  More  specifically  we  wish  to  show  that 


P  (x) 
n 


9 


uniformly  in  every  circle  of  IT. 

In  1913  F.  Riesz  devoted  a  hook  [6]  to  such  problems  which  he  calls 
"Problfemes  des  redultes".  His  general  theory  does  not  apply  to  our  specific 
problem.  However,  further  work  in  Functional  Analysis  might  solve  it.  This 
shows  that  Fulerian  mathematics  still  presents  to  contemporary  analysts 


challenging  problems. 
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n 
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